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Abstract 

J/iJj and ip 1 production cross sections in fixed-target experiment is calcu- 
lated, considering the contributions from both color-singlet and color-octet 
mechanisms. The results are applied to the investigations of the J/ip sup- 
pression and the ip' /ip ratio problems in p-A collisions. The results agree 
with the experimental data as the (cc)— nucleon absorption cross sections 
a^ bs ~ Wmb for (cc)s and a\ bs ~ Omb for (cc)i. The model is further used 
to investigate A-A collisions when comover absorption mechanism is also 
considered. It is found that the observed experiment data of J ftp and ip' /ip 
ratio in S-U collisions and Pb-Pb collisions can not be explained consis- 
tently within this model. The possibility of QGP formation in S-U and 
Pb-Pb collisions is also discussed. 
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I. Introduction 

Matsui and Satz proposed that a suppression of J /if) production in relativistic heavy 
ion collisions can serve as a clear signature for the formation of a new matter phase Quark 
Gluon Plasma (QGP) |T|. This suppression effect was observed by NA38 collaboration 
later |2[]. However, it has been found that J/if> suppression exists also in p-A collisions 
where QGP formation is not possible ||. The successive theoretical researches pointed out 
that J /if) suppression could also exist in hadronic matter (HM), even though by completely 
different mechanisms ||. The anomalous J /if) suppression was recently reported by the 
NA50 collaboration || || and there have been a number of attempts to explain it, such 
as the onset of deconfinement, hadronic co-mover absorption and the energy loss model 
@ || 0. To understand the experimental data clearly, the formation and absorption 
mechanism of J/ if) must be studied carefully. 

Quarkonium production has traditionally been calculated in the color singlet model. 
However, it has become clear now that the color singlet model fails to provide a theoretical 
and phenomenological explanation of all J /if) production processes and it is necessary to 
include the color octet production channel. In principle, the J /if) state is described in a 
Fock state decomposition 

\J/$ > = 0(l)\cc( 3 S{ 1] ) > +0(v)\cc( 3 Pf)g > +0(v 2 )\cc( 1 S^ ) )g > + 

0(^)\c^S^)gg > +0(v 2 )\cc( 3 D { j'%g > +..., (1) 

where 2S + 1 £ l y 8 > characterizes the quantum state of the cc in color-singlet or octet |10| |, 
respectively. This expression is valid for the non-relativistic QCD (NRQCD) framework 
and the coefficient of each component depends on the relative three- velocity \v\ of the 
heavy quark. Under the limit of \v\ — > 0, i.e. c and c remain relatively at rest, Eq.(|l]) 
recovers the expression for color- singlet picture of J /if), where 0(1) = 1. 

The color-octet component (cc)§ seems to play an important role in interpreting the 
experimental data by Collider Detector at Fermilab (CDF ) jnj, and further studies show 
that it may also have great influence on the quarkonium production at other collider facili- 



ties [pTT[— 113[| . The investigations |14| on quarkonium hadrpro duct ion at fixed target energies 
find that the color-octet contribution to the production cross section is very important, 
and the inclusion of color-octet production channels removes the large discrepancies of 
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the total production cross section between experimental data and the predictions of the 
color-singlet model. Therefore, it has been suggested that the color-octet (cc) 8 would also 



manifest itself in heavy ion collisions f[? |. Based on the above discussion, one can accept 
the following physical picture that the charmonium production can be divided into two 
steps. The first step is the production of a cc pair. The cc pairs can be either (cc)i or 
(cc)s, which are produced perturbatively and almost instantaneously, with a formation 
time t f ~ (2m c ) _1 ~ 0.07/m in the cc rest frame. The second step is the formation of 
a physical states of J/ip, that needs a much longer time. People believe now that J/ip 
suppression in hadron matter can be considered as the pre-resonance absorption. Satz 
first proposed the pre-resonance absorption model to explain the nuclear collision data 
and got some encouraging results [|17]]. However, in their work the pre-resonance state of 
charmonia is only in color-octet, the color-singlet counterpart is believed having a minor 
influence, therefore negligible. In this paper we make a complete leading order calculation 
of charmonium production cross section for (cc)i and (cc)s in concerned energies. The 
results are used to study the charmonium suppression in heavy ion collisions. 

As it is known, in hadron-hadron collisions the initial (cc) pair produced from the 
parton interaction are almost point like. It can be in color-singlet or -octet configura- 
tions. The (cc) pair in the color-octet state may interact with the nucleon environment 
much stronger than that in the color-singlet, which means that the color-octet pair would 
be dissolved in a much shorter time than the color-singlet pair. Based upon this, the 
absorptions of (cc)i and (cc)s by nucleon environment should be considered differently. In 
the framework of this paper, we suppose that the (cc)— nucleon absorption cross sections 
al bs for (cc)i and a^ bs for (cc)s satisfy the condition a\ bs ~ << <rf bs , however the color- 
singlet part is not completely negligible. With this assumption the J/ip suppression data 
in p-A collision can be explained. The experimental data of tp' /if) ratio in p-A collisions 
can also be explained quite well. For nucleus— nucleus (A-A) collisions, in addition, the 
comover absorption is investigated as well and we find that the experimental data of J / ip 
suppression and ip' /ip ratio can not be described consistently in S-U and Pb-Pb collisions. 
The possibility for the production of QGP is also discussed. 

Our paper is divided into four parts. In section II, we describe how to introduce the 
color-octet scenario borrowed from p-p collisions and calculate the charmonium produc- 
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tion cross section through both color-singlet and color-octet mechanisms. In section III, 
the pre-resonance nuclear absorption model relating to both (cc)i and (cc)§ is described 
and used to discuss the J/ip and ip' suppression in p-A collisions. In section IV, the 
nucleon and comover absorption model for J/ip suppression is investigated. In the last 
section, some further discussions are given. 

II. Formulation 

According to the NRQCD factorization formalism, the inclusive production rate of 
heavy quarkonium H in parton level can be factorized as 

a(ij -#) = £ - QQ[n]) < H [n] > . (2) 

n 

Here, a(QQ[n] +X) describes the short distance production of a QQ pair in the color, spin 
and angular momentum state n, which can be calculated perturbatively using Feynman 
diagram methods. < (9 H [n] >, the vacuum expectation value of a four fermion operator 
in NRQCD |1(J, describes the nonperturbative transition of the QQ pair hadronizing into 



the quarkonium state H. The relative importance of the various contributions of n in 
Eq.(0) can be estimated by using NRQCD velocity scaling rules. An important feature 
of this equation is that QQ pairs in a color-octet state are allowed to contribute to the 
production of a color singlet quarkonium state H via nonperturbative emission of soft 
gluons. Accordingly, the production cross section for a quarkonium state H in the hadron 
process 

A + B -> H + X (3) 

can be written as 

<*H = £ / dx 1 dx 2 fi/A(x 1 )f j /B(x 2 )cr(ij -> H), (4) 

where the parton scattering cross section is convoluted with parton distribution functions 
fi/A and fj/Bi and the sum runs over all partons in the colliding hadrons. 

At leading twist and at leading order in a s , the color-singlet QQ production subpro- 
cesses for 2S+1 Lj state are 
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gg - 3 P , 2 , (5) 

3 Si+£, 3 Pi + £, (6) 

<?g^ 3 P + g, (7) 

gg^ 3 Pi+#. (g) 



The corresponding formulae of the above precesses can be found in refs. |14| and [15 ]. 
Because the radiative decays xi,2 ~^ J/ip + 1 are known to have a large branching ratios 
to J/ip and the feeddown of the ip' to J ftp is also important, their contributions should 
be included in the calculation in reproducing the fixed target experimental data of the 
prompt J/ip production. 

It has long been known that the total cross sections of quarkonum production are 
rather large in comparison with the fixed-target experiment data with respect to the 



parton-parton fusion predictions JTB . However, after including the color-octet production 



mechanism an overall agreement can be obtained |14 



Under the NRQCD factorization scheme, to calculate the quarkonium production cross 
section, one use a double expansion: the perturbative expansion of the short distance 
production amplitude in strong coupling constant a s and the expansion of the nonpertur- 
bative long distance hadronization amplitude in typical velocity of heavy quark inside the 
heavy quarkonium. At leading order in perturbative theory and up to next-to-leading or- 
der in the velocity expansion, the subprocesses for leading-twist J/ ip production through 
color-octet intermediate states are 

qq -> cc[8, 3 S 1 ] -> J/V> + X, (9) 
gg^cc[8, 1 S ]^J/ij + X, (10) 
^^cc[8, 3 Pj] - J/^ + X, (11) 
qq -> cc[8, 3 Pj] -> XJ + X - J/i/j + 7 + X. (12) 

Note that the process qq — > cc[8, 3 Pj] — > J/ip{ip') + g is of higher order in v 2 since 
the lowest-order nonperturbative transition if forbidden by charge conjugation, and the 
amplitude A(gg — ► ccfS, 3 ^]) vanishes in leading order in at s , which is consistent with 
Yang's theorem which forbids a massive J = 1 vector boson from decaying to two massless 
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J = 1 bosons In fact, the theorem requires that the A(gg — > cc[8, 3 Si]) vanish to all 
orders as the gluons on the massshell. 

The cross sections of Eq.(||) are proportional to the NRQCD matrix elements 

<O|08^( 3 Si)|O>~ rn 3 c v 7 , (13) 

<0|O 8 J/,/, ( 1 5 )|0>~my, (14) 

<0|O^( 3 Pj)|0>~my, (15) 

< 0\C$ J ( 3 Si)|0 >~m 3 c v 5 . (16) 

It is obvious that the above matrix elements are higher order in v 2 compared to the 
leading color-singlet ones, but their corresponding short-distance processes are lower order 
in a s than that in color-singlet processes. This causes the color-octet processes to make 
substantially enhancements in reproducing the fixed target experiment data. 

For ip' production the cross section does not receive contributions from radiative de- 
cays of higher charmonium states. <j(ip') differs from the direct J ftp production cross 
section u{J/ip)dir only in the replacement of J/if> matrix elements in Eq.fll3|) by ip' matrix 
elements. 

Before embarking on the computation of cross sections the parameters used in the com- 
putation should be fixed up. The uncertainties in the theoretical prediction at fixed-target 
energies are substantial and it is impossible at present to extract the non-perturbative uni- 
versal color-octet matrix elements by fitting the theoretical predictions to the experiment 
data. The value of < 0|Cg^^ ^( 3 Si)|0 > used here is obtained by fitting the theoretical 
predictions to the CDF Collaboration data at large pt- In addition, the number of in- 
dependent matrix elements can be reduced by using the spin symmetry relations up to 
corrections of order v 2 

< O|0? J ( 3 Pj)|O >= (2J+ 1) < O|0f ( 3 P )|0 >, (17) 
< 0|O^( 3 Pj)|0 >= (2J+ 1) < 0|O 8 W ( 3 P )|0 >, (18) 
<O|0£ J ( 3 Si)|O>= (2J+1) <0|Of( 3 5'i)|0> . (19) 

Therefore, the matrix elements < 0|C , 8 ¥ ( 1 S , )|0 > and < 0|O 8 ¥ ( 1 P )|0 > enter fixed target 
production of J/ip and ip' in the combination 
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A 8 (H) =< 0|Of (^0)10 > +\ < O|0f ( 3 P )|0 > . (20) 

Tfls 



Up to corrections in v 2 , the relevant color-singlet production matrix elements are related 
to radial wave functions at the origin or their derivatives, 

< 0|Of ( 3 ^)|0 >= ^|P(0)| 2 , < O|0f ( 3 P )|0 >= ^-|i?'(0)| 2 , (21) 

which can be determined from potential model or from quarkonium leptonic decays. 
The values of these parameters, which we use, are [[L4j] 

< O|0J w ( 3 S , i)|O >= 1.16 GeV 3 , < Opt'* ( 3 Si)\0 >= 6.6 x 10" 3 GeV 3 , (22) 
< 0\Of ( 3 5i)|0 >=0.76 GeV 3 , < 0\Of ( 3 Si)|0 >= 4.6 x 1(T 3 GeV 3 , (23) 

< 0|Of ( 3 P )|0 > /m 2 c = 4.4 x lO- 2 GeV 3 , < O|0£°( 3 Si)|O >= 3.2 x 1(T 3 GeV 3 , (24) 

A 8 (J/iP) = 3.0 x 10" 2 GeV 3 , A 8 (^') = 5.2 x 10" 3 GeV 3 . (25) 

In the numerical calculation, we use the Gluck-Reya- Vogt (GRV) leading order (LO ) 
|20| parameterization for the parton distributions of the protons. The c quark mass is 
fixed to be m c = 1.5 GeV and the strong coupling is evaluated at the scale \l = 2m c , 
that is a s ~ 0.26. The results of the integrated cross sections for color singlet and color 
octet channels at several different energies are listed in Table I. We must admit that the 
results are far from precision because of some important contributions may be precluded 



as claimed by authors [14| [21|. e.g., the effects of higher twist, the beyond leading order 



contributions in a s , as well as the kinematic, etc.. 

III. Pre-resonance Absorption in p-A Collisions 

The conventional survival probability for a J/ip produced in a p-A collision is given 

by: 

1 CTpA 



S A 



A dpp 



d bdzpAip, z)exp< —(A — 1) / dz' ' Pa(P, z')a, 



' abs 



exp{-L A po<J abs ), (26) 
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where a pp and <j v a are the J /if) production cross section in proton-proton collisions and 
proton-nucleus collisions, respectively, pA is the nuclear density distribution. cr a b s is the 
absorption cross section. L is the effective length of the J/ip trajectory. It can be derived 

as 

3 A — 1 

L = — — — — tqA 1 / 3 , for heavy nucleus 
1 A — 1 i r 2 

= — — — — roA 1 ' 3 -^, for light nucleus, (27) 

where po — 0.14/m -3 and ro = 1.2/m, r f = 1.05/m. 

As cc pairs are produced almost instantaneously and the formation of the physical 
states J/ip or ip' need a much longer time, people now believe that J/ip and ip' suppression 
in p-A can be considered as an absorption of pre-resonance cc pairs. As discussed in former 
section, there are both (cc)i and (cc)s pairs. The color-octet can interact with gluons much 
more strongly than the color-singlet (cc)i, and therefore would dissolve much faster into 
D and D than (cc)x- Thus their absorption cross section are different. Considering these 
facts, we rewrite Eq.(^) as 

1 &pA 



A o~p P 



Sa = 

= h J d 2 bdzp A (b,z)exp^-(A-l) J dz' p A (b, z')al hs 

+ f 8 J d 2 bdzp A (b, z)exp{-(A - 1) dz'p A (b, (28) 

where /i,/s are relative fractions of (cc)i and (cc)s- c^lbs^abs are ^ ne absorption cross 
sections for (cc) i-nucleon and (cc) g-nucleon, correspondingly. 

In Eq.(|2"%D there are two parameters, a\ bs and a^ bs , which is different from Satz's model 



T?]]. As (cc) 1 produced is almost point like, a\ bs is very small. Thus we take a\ bs = and 
the value of a® bs is considered as an open parameter and determined such as to get the 
best agreement with the data. In Fig.l we see that with o~l bs = lOmb we get quite good 
agreement with the p-A data. 

Next we turn to discuss the ratio ip'/ip in p-A collisions. As J/ip and ip' suppression 
in p-A can be considered as an absorption of pre-resonance cc pairs, there is no difference 
for J/ip and ip 1 in p-A collisions. Then the ip ! /ip ratio in p-A collisions can be expressed 
as 
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B{J/lfj -> p+p^Vp-A^J/V 



(29) 



B(J/V> -> //+// )[(Tiexp(-L A p ^D + ff 8exp(-L A p ^D] ' 
where B is the corresponding branch ratio, a±, a 8 are the production cross section of color 
singlet and color octet for J /if) in p-p collisions, respectively. a[, <r' 8 are the production 
cross section of color singlet and color octet for ip' in p-p collisions, respectively. 

Using the same parameters as those in Fig. 1, the result of Eq. ( ^9]) is shown in Fig. 2, 
where one can see that the ratios ip' /ip obtained in our model, being almost independent 
of the c.m.s. energy, agree with the experimental data quite well. The results show that 
including both the contribution of color singlet and color octet does not introduce any 
unusual A dependence. 

IV. Comover Absorption in A-A Collisions 

In A-A collisions, except for (cc)— nucleon absorption, charmonium may also suffer 
interaction with secondaries that happen to travel along with them. The J /if) survival 
probability due to absorption with comover hadrons is 

S°° = exp{- J dra co n co (T, b)}, (30) 

where a co is the J/^-comover absorption cross section, n co (r,b) is the density of the 
comovers at time r and impact parameter b. The relative velocity between J/if) and the 
comover is included in the definition of the absorption cross section. 

Integrating over time r, assuming that the comovers undergo an isentropic longitudinal 
expansion, Eq. fl3"0|) can be expressed as 

S co = exp{-a co n T ln—}, (31) 



S m = exp{-a co n T ln—}, (32) 

where tq is the production time of the comovers and t$ is the time the comovers and 
J /if) stay together, no and rif are the initial density and freezout density of the comovers 



separately. When is smaller than the life time of the comovers t^, Eq.(^) describes the 
comover survival probability, otherwise Eq.(|3l|) works. For if)', there are similar equations 
except that o co is replaced by a' co which is the -^'-comover absorption cross section. As 
the mass of if)' is much closer to the DD threshold, only a 50MeV excitation is needed to 
break up a if)', while for J /if), nearly 650MeV is needed to be above the DD threshold. 
Thus a' co should be much larger than a co . 

Considering the effects of both the (cc) — nucleon absorption and J/^-comover absorp- 
tion, the J /?p survival probability in A-B collisions is 

S = S co x S nuc , (33) 

where S nuc is the J /if) survival probability in A-B collision due to the (cc)— nucleon ab- 
sorption. It is similar with Eq. (p8|) and can be expressed as 

S nuc = f x exp(-(L A + L B )p al bs ) + f s exp(-(L A + L B )p a 8 abs ), (34) 

where La and Lb are the effective length of the J /if) trajectory in A and B nucleus 
correspondingly. 

From the above equations, one can see that there are some parameters in the comover 
model. In Eq . (|3l"D , the parameters are: n , a co , r , rij. In Eq . (|32|) , the parameters are: n , 



Ceo, tq,^. We first discuss the data at different Er-bins and analyze the Et dependence of 
the parameters. In this paper we wish to adjust these parameters consistently for S-U and 
Pb-Pb collisions. The comover density is taken to be proportional to the energy density 
at a certain space-time point, which can be expressed as the density of the transverse 
energy according to Bjorken's assumption ||22|| . Therefore, the comover density could be 
expressed as 

E T E t 

"° ~ aF ~ swaw (35) 

where AV is the corresponding volume. S(b) is the overlapping area of the two nuclei. 
Ay is the corresponding rapidity windows in the central rapidity region. Ay = 2.4 for 
S-U collisions and Ay = 1.2 for Pb-Pb collisions. 

Considering Eq.(|35D the Et dependence of the comover density in S-U and Pb-Pb 
collisions can be described at the same time based on collision geometry. Now we use 
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Eq.([33l) and Eq.([3~Tl) to fit the experimental data in S-U and Pb-Pb collisions. The pa- 
rameters a\ bs and a^ bs used in considering the (cc)-nucleon absorption in A-B collisions is 
taken to be the same as those obtained in fitting the data in p-A collisions. We treat the 
comover density in the first Et bin of S-U collisions, nj, as an open parameter, adjust 
nl,nf,cr co to fit the J /if) suppression data in S-U collision, then choose a' co to get the 
best fit for if>'/if> ratio data in S-U and Pb-Pb collisions. The results for S-U and Pb-Pb 
collisions are shown in Fig.3-4, with parameters = 0.2/m~ 3 , o" co = 3mb,a' co = 23m6, 
and nf = 0.1/m -3 . 

If t$ <t^p, one should use Eq. (|32"D to describe the comover absorption. It is reasonable 
to choose T^, proportional to the square root of the transverse overlapping area of the two 
nuclei 



~ ft®. (36) 



Using Eqs. (|35|) , (j36|) , (|32|) and (|33|), with the parameters = 0.2/m~ 3 , a co = 3mb, a' co = 
I3mb, and ri = 6.5/m, which is the for the first Et bin in S-U collisions, the results 
are shown in Fig. 5 and Fig. 6. 

Fig. 3-6 show that neither of the two comover absorption expressions of Eq.(|31]) and 
Eq.(|52"D can explain the data of if)' and J /if) suppression in S-U and Pb-Pb collisions 
consistently. Fig. 3(a) shows that the J /if) suppression in S-U collisions could be described 
very well, based on the above chosen parameters for Eq.(|3*ID, however Fig.3(b) shows that 
the if)'/ if) ratio data in S-U collisions could not be fitted using the same set of parameters. 
The data show an anomalous if)' suppression from the second Et bin. Furthermore, 
Fig. 4(a) shows that the comover absorption which explains the J /if) suppression in S-U 
collisions can not explain the J /if) suppression in Pb-Pb collision, where an anomalous 
suppression exists from the second Et bin. This seems to show that the anomalous 
if)' suppression begins already in S-U collisions at an energy density much lower than 
the corresponding density for anomalous J /if) suppression in Pb-Pb collisions. This may 
reflect the fact that the dissociation temperature for J /if) in QGP is higher than the 



dissociation temperature of if)' | 23fl . From Fig.4(b) one can find that the same comover 
absorption can explain the data of if>'/if) ratio in Pb-Pb collisions. This could be explained 
as that the anomalous J /if) and if' suppressions are canceled in the if' /if ratio data. The 
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results of using Eq. (|32"D to include the comover absorption are shown in Fig. 5-6, which 
are similar to Fig. 3-4. 

Now we turn to discuss the case of minimum biased data, where the result is 
shown in a simple and clear way. With the parameters rif — 0.1/m -3 , tq = lfm, 
o co = 3mb, o' co = 23mb which is the same as that we used in obtaining Fig. 3-4, taking 
the average comover density n su = 0.28 in S-U collision and the average comover density 
n p b = 0.4 in Pb-Pb collisions, the results of using Eq.(|33"D and Eq.([H]) f° r the minimum 
biased data are shown in Fig. 7. From it one can see clearly that using the comover 
absorption expressions of Eq.(|3lD the J /if) suppression data for S-U collisions is fitted 
based on the above parameters, but one can not explain the data of if)' suppression in S-U 
collisions, neither the J/ if) suppression in Pb-Pb collision. The good fitting of if}'/ if) ratio 
data in Pb-Pb collisions may caused by the same anomalous absorption of if>' and J/if>. 
The results of using Eq. (|33|) and Eq.(32) is similar to Fig. 7. 



V. Results and Discussions 

In this paper, J/if> and if>' production cross section is calculated considering the con- 
tributions of both color-singlet and color-octet cc channels. The pre-resonance absorption 
model for charmonium is extended to consider both the color singlet and color octet con- 
tribution. Using this model the J/ if) and if}' suppression in p-A collision are explained very 
well. Based on above calculation, the comover absorption is discussed for A-A collisions 
and it is found that the observed experimental data of J / if> and if)' / if) ratio in S-U collision 
and Pb-Pb collision can not be explained consistently by this mechanism. This indicates 
that other sources of charmonium suppression should be included. 

The situation in explaining the data of strangeness production is similar. To explain 
the enhanced production of strangeness, some kind of collective interaction among indi- 
vidual excited nucleon states, such as the colour ropes in RQMD |24j, the multiquark 



clusters in VENUS ||25|| , firecracker in LUCIAE p6| , must be considered. Some authors 



also have reported the possibility of QGP formation at CERN SPS from the strangeness 
puzzle |27j . There is no strict way to distinguish these different kinds of collective motions 
from the formation of QGP till now. 
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Now we consider the possibility of QGP formation from our discussion. In QGP, 
charmonium breaks up due to color screen. As the radius of ip' is larger than that of J/ip, 
the critical temperature T' c for ip' to be dissociated in QGP should be much lower than 
the critical temperature T c for J/ip to break up [[^J. So if QGP is formed, if)' will begin to 
break up earlier than ip. As has been pointed out above, Fig. 3(b) and Fig.5(b) indicate 
that from the second Et bin (or the third Et bin) in S-U collisions there is an anomalous 
ip' suppression, which may be the effect of QGP production. Fig. 4 (a) and Fig. 6 (a) show 
that the anomalous suppression of J/ip really begins after the second Et bin in Pb-Pb 
collision because it needs a higher critical temperature. If deconfmed phase is attained, 
when J I ip begin to be suppressed due to Debye screen, ip' should already be broken up by 
Debye screen. These two anomalous suppressions in QGP may cancel each other in ip' /ip 
ratio, which leads to the result that ip'/ip ratio in Pb-Pb collisions seems to be explained 
by the comover absorption, while in fact, neither J/ip nor ip' anomalous suppression could 
be explained by comover absorption. Fig.4(b) and Fig. 6(b) agree with the above picture 
of J/ijj and ip' suppression by QGP. 

Although our results indicate the possibility of production of QGP in S-U and Pb-Pb 
collisions, and to explain enhanced strangeness production some collective motion which 
can not be distinguished from QGP formation must also be considered, further study is 
still needed to distinguish if QGP are really formed. 
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Figure Captions 

Fig.l: The J/ip survival probability obtained using Eq.(|2"%|) with a^ bs = lOmb and 
a lbs = * s compared to the experimental data of p-A collisions at different energies. 

Fig. 2: The ip'/ij) ratio obtained using Eq. with a^ bs = lOmb and cr\ bs = Omb is 
compared to the experimental data of p-A collisions at different energies. 

Fig. 3: (a) J/ip over DY ratio versus E T and (b) B ^cr^') / B ^(^(b) versus E T in 
S-U collision are compared to the results based on Eq . ([33D using the comover absorption 
Eq.(pT|). The parameters are n\ = 0.2/m~ 3 , a co = 3mb, a' co = 23mb, and rif — 0.1/m -3 . 

Fig. 4: The same as Fig.3 for Pb-Pb collisions. 

Fig. 5: The same as Fig.3 using the comover absorption Eq.(p2[) with the parameters 
n\ = 0.2/m -3 , a co = 3mb, a' co = 13mb and ri = 6.5 fm. 

Fig. 6: The same as Fig.5 for Pb-Pb collisions. 

Fig. 7: The same as Fig. 3 in the minimum bias case with n su = 0.28, n p b = 0.4. 



Table Caption 

Table I. The integrated cross sections for color singlet and color octet processes. 
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Table I 
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